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If the universe is trapped and cooled in a metastable false vacuum state, that state will eventually 
decay by bubble nucleation and expansion. For example, many extensions of the standard model 
incorporate new scalar fields (f> whose potential has a local minimum at <f> = but a global minimum 
elsewhere, to which the vacuum will eventually tunnel. I calculate the lifetime of the false vacuum, 
and the field profile of the bubble after tunneling, for any potential that is approximately a polyno- 
mial of degree < 4 near the false vacuum. Essentially exact results are given for a single field; for 
multiple fields a strict lower bound is placed on the tunneling rate. 
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Consider a quantum field theory of some scalar fields <p 
(grouped for convenience into a vector) whose potential 
U(<f>) has local minima at the origin <f> = and at some 
other (j> = <f>\. If U(<pi) < U(0) then the vacuum at the 
origin is a false and metastable one. If the universe is in 
this false vacuum at some early epoch, and the tempera- 
ture is sufficiently low - as I will assume for the remainder 
of this paper - then quantum fluctuations will eventually 
initiate a phase transition from the supercooled initial 
state to the true vacuum through bubble nucleation and 
expansion [EJ. A very elegant and tractable approach 
to studying this process was introduced by Coleman [Q. 
The bubble's most likely field configuration, called the 
bounce, is the one which extremizes the Euclidean ac- 
tion, and the lifetime of the false vacuum is proportional 
to the exponential of this extremal action. The bounce 
also determines the bubble's form in Minkowski space, 
after its formation and outside the lightcone of its cen- 
ter. 

In this work I calculate the bounce and its Euclidean 
action for a wide class of models, namely all those for 
which the potential is at least approximately a polyno- 
mial of degree < 4 near the false vacuum. The results 
are essentially exact for a single field <p; with more fields 
the action and therefore the false-vacuum lifetime I cal- 
culate are strict upper bounds, while the one-dimensional 
bounce given here may only give a qualitative picture of 
the true field configuration. One class of applications for 
this work is to extensions of the standard model contain- 
ing fields <f> whose vacuum expectation values (VEVs) 
in the present epoch are phenomenologically required to 
vanish, for example if they are electrically charged. Thus 
if a lower minimum of U develops away from the origin, 
the origin must not only remain a local minimum, but 
must also be sufficiently long-lived that we would prob- 
ably still be inhabiting it today. For this type of appli- 
cation the lifetime is of primary importance. It is hoped 
that, within other contexts, the shape of the bounce will 
be useful in calculating various other properties of the de- 
cay, for example the degree of supercooling, the thickness 
of the transition region and the distribution of energies 



within its volume. Of course, for many of these con- 
texts the effects of gravity H and the expansion of the 
universe, and possibly of finite temperature B, must be 
included, but the techniques developed here may readily 
be extended. 

Most previous work on false vacuum decay either as- 
sumes a thin-wall approximation to derive analytical re- 
sults, or specializes to a particular model and numerically 
extremizes the action. But the former approach is often 
only applicable when the lifetime is extremely long and 
hence not very interesting, while the latter can be intri- 
cate and time-consuming, requiring a new computation 
for each model and often also ingenious methods of ex- 
tremization ||. The results presented here yield either 
a bound on the action or an essentially exact answer for 
a very wide class of potentials, with any choice of pa- 
rameters and without any computer reanalysis. There 
is nevertheless some overlap with previous authors. For 
the case of a single scalar field cf>, rescaling was used (as 
I do below) to exactly study two limiting-case potentials 
HQ and a general quartic potential though most 

of the latter results are given only graphically and for 
a limited range of parameters; in the high-temperature 
limit the action was calculated and presented completely 
in Ref. [Eo|. For the case of multiple fields $, the reduc- 
tion to a single field was employed, for example, in Refs. 
||s|,fi"l| . These were particular realizations of some of the 
same techniques used in the present work, and, as further 
emphasized in Ref. fl^| , illustrate the usefulness of this 
approach. 

Consider any potential U in any number of scalar fields 
<f> with a homogeneous false vacuum at <f> = and with 
U(0) = 0. In the semiclassical approximation 0, the 
probability of bubble formation per unit time and volume 
is T/V ~ to 4 exp(— Se)- The prefactor is dimensionally 
the fourth power of the typical mass scale in the theory 
and thus can be readily estimated. I will concentrate ex- 
clusively on the argument of the exponential, which is 
the Euclidean action for the bounce solution. However, 
care should be taken to account for any large numerical 
prefactors, for instance p3 group-theoretical factors re- 
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lated to vacuum degeneracy in models with spontaneous 
symmetry breaking. 

The requirement that the false vacuum has probably 
not decayed in our past light-cone reads: (T/V)L 4 <C 1, 
where L is roughly the present size and age of the vis- 
ible universe. Numerically, L ~ 10 10 yr ~ 10 45 TeV _1 
so S E > 400 + 41n(m/TeV). The bounce is an 0(4)- 
symmetric JT3] solution of the Euclidean equations of 
motion which approaches the false-vacuum field configu- 
ration <j) — at Euclidean time r — ±oo. To calculate 
it, one must extremize the action with respect to all pos- 
sible non-trivial paths through field space which satisfy 
the boundary conditions. This requires some careful and 
somewhat time-intensive numerical techniques and 
the results are difficult to generalize. To simplify this 
task, consider a single, straight-line path: cf> = V(f> where 
v is a constant. For example, v can be fixed to point 
towards the true vacuum, or towards any other direction 
where tunneling is possible. (The action can also later 
be minimized with respect to v, or to small perturba- 
tions thereon [0). In any case, the extremum of the 
true action Se must be no larger than the extremum of 
the action S% fo r (f> restricted to lie along v, by a theorem 
of Coleman [[15| . Any model not obeying 



S* E > 400 + 41n(m/TeV) 



(1) 



will have an even smaller Se and hence its false vac- 
uum would have decayed by today. Constraining a model 
based on Eq. ([!]) isn't as stringent as contraining it based 
on the full Se, but it is certainly much simpler, can be 
studied in general as I do below, and is often sufficient to 
rule out large regions of parameter space. For instance, 
Dasgupta jl2j has shown for a particular supersymmetric 
model that most of the parameter range excluded by an 
exhaustive Se extremization is already ruled out by the 
simple Se constraint; see also Ref. ||. 

Thus we consider a potential function of a single field 
<f>, in particular: 



U{(j>) = Aff 2 - M 3 (j) 3 + \(f) 4 



(2) 



This form encompasses all renormalizable potentials, but 
is also a good approximation for many nonrenormalizablc 
ones which can be expanded as a quartic polynomial near 
the false vacuum. I assume M| > to ensure at least 
metastability of the origin, and M3 > without loss of 
generality. By Euclidean spherical symmetry, the action 
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FIG. 1. The rescaled Euclidean action Se as a function of 
k, with asymptotes as dashed lines. 



Se = 27r" / ,c" 
'0 



dx 



U(y) = y 2 -y 3 + Ky 4 , k = AAff /Mf . 



(4) 



(5) 



Thus the general problem has been reduced to extrem- 
izing an action Se depending on a single parameter k. 
Tunneling is possible for any | > k > —00; when n < 
some other mechanism must eventually stabilize the po- 
tential, but will not influence the following if the true 
potential is well- approximated by U at all < y < yo, 
where yo is the escape point defined below. 

Extremizing Se amounts j^j to solving the equations 
of motion 



d 2 y 



dy 



| 3dy _ d(-U) _ 
dx 2 x dx dy dx 



= 0, 2/(00) = (6) 



for a particle located at position y as a function of time 
x > and subject to a potential — U[y(x)] and a time- 
dependent viscous frictional force. The particle starts at 
some y{x = 0) = j/o > and slides down the potential 
— U monotonically towards the origin, with yo chosen so 
y(x — > 00) = 0. 

The above equations are numerically integrated 
throughout the allowed range for n. The result- 
ing action, from Eq. (^) or equivalently from Se — 
— 2tt 2 J °° x 3 U[y(x)] dx, is shown in Fig. [j]. To better than 
1% accuracy this action is fit by the semiempirical expres- 
sions 



S E = 27i- 



zd^4> + u{4>) 



dp 



(3) 



where z is some constant factor resulting from field rescal- 
ing (z = 1 is canonical normalization) and p = Vx 2 + t 2 . 
Using a dimensionless coordinate x = (M2/ z 1 ' 2 )p and 
field y(x) - {M 3 /M 2 ) <f>(p), I find S E = (z 2 M 2 /M 2 ) S E , 
where 



Se,+ — Se, thick — 46.1 + S 



B.thin 



16.5 28 

+ - r- (7) 



Se,- — Se, thick 1 + (Se, thick/ Se,f) 1a N 1 



(1-4k) 2 I -4k 



(8) 



where the + (— ) subscript indicates k > (k < 0), 
Schick =s 45.4 (see also §), S E ,tUn = 2tt 2 /[12(1 - An) 3 } 
and Se.f = 27r 2 /3 ensure the correct action in the thick 
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FIG. 2. The (rescaled) bounce y(x) for 0.05 < x < 20 
and various k values: (A) -5, (B) -2, (C) -1, (D) -0.5, (E) 
-0.25, (F) -0.1, (G) -0.01, (H) 0.01, (I) 0.1, (J) 0.175, (K) 
0.2, (L) 0.22, (M) 0.23, (N) 0.235. 
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FIG. 3. The (rescaled) bounce is well-approximated by 
VRt(x) = yo/[l + x 1 e^ x ~ R) ]\ the "initial" value y and the 
fit parameters R and 7 are plotted in (a), (b), and (c), re- 
spectively. 



wall, thin wall and Fubini asymptotic limits, respectively, 
as discussed below. 

The bounce solutions to the equations of motion are 
plotted in Fig. || for a various k. Perhaps surprisingly, 
they all fit quite accurately the simple form 



y(x)~y Rt (x)=y / [l + ^e^"^ 



(9) 



where y = y(x = 0), while 7 and R are empirically 
chosen as functions of k to optimize the fit. Fig. || shows 
yo and the best-fit parameters 7 and R as functions of k. 
Note that when k is positive but significantly below its 
maximal value, say < 0.15, the shape of the bubble and 
especially its action are very different from the thin-wall 
results. In particular, the field tunnels out to a value yo 
far from the true vacuum value §k _1 [l + (1 — if k) 1 / 2 ] 
shown as the long-dashed gray curve in Fig. ||a. The fit 
parameters are themselves well-approximated by 



yo,+ =s ?yo,thin + 1-18(1 - 4k) 2 - 1.40(1 - 4k) 3 
yo,- ~ yo,thick + 10.4 ln(l + \k\) - 1.67 ln(l + k 2 ) 

7+ ~ 2- [2- 7.2k- 11.4k 2 - 0.72 ln(l - 4k)]" 1 

7 _ ~ 7f -1/[2 + 4.9|k|-1.3|k| 1 - 3 ] 

R + ~ R thin + 0.91/(1 - 4k) 1/2 - 1.41 + 6.4k 12 

J?_ ~ 0.21 - 1.13 ln(l + 11|k| 115 ) 



(10) 

(11) 
(12) 

(13) 
(14) 
(15) 



Here y ,thin = 2 + 4(1 - 4k), y ,thick ^ 5.78, 7f = 2, 
and iZthin = [v2(l — 4k)] _1 are the asymptotic values 
described below. The y/2x in the exponential makes 
2/fit tend towards zero at the correct rate as x — > 00, 
and the various expressions (|Tc|)-(15) give an excellent 
fit wherever y(x) is significantly different from zero: 
\ytit{x) - y{x)\/y ^ 1% for all x. 

Several limiting cases are of interest. The best-known 
is the thin- wall regime |Q, when k — > 1/4 and the 
two vacua approach degeneracy. To leading order in 
1 — 4k the equations of motion may be solved analyt- 
ically: ythin = 2/o,thin/{l + exp[ v / 2(x - i?thin)]}- Note 
that, while in general the solution is better fit by Eq. (|^) 
with finite 7, in the thin-wall limit 7 becomes entirely 
immaterial and ym — > ythin- The opposite limit, in a 
sense, is when k — ► 0, which I call (as did Ref. [fl6]) 
the thick- wall regime. The values 7thick — 1-5 and 
-Rthick — 0.21 provide a better fit over the entire range 
of x than the exact solution of Ref. jl6) to the linearized 
equations of motion. A more subtle limit is k — * —00, 
discussed in some detail in Refs. I^Q. Then the cu- 
bic term in U becomes irrelevant, so Se — * \k\~ 1 Se,f 
where Se,f = 27r 2 J °° x z \\(dyF / dx) 2 + y 2 F — y F ]dx and 
jjf = l^^y. This action can only be extremized asymp- 
totically. Consider the action of the family of functions 
y Ff} = V2/3/(/3 2 + x 2 ). When /3 -> 0, the quadratic 
terms may be neglected, in which case the yF\f3 are ex_ 
act solutions (known as Fubini instantons |17| ) with a 
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/?-independent action |-7r 2 . Thus when the complete ac- 
tion Se,f is computed for yp,/3 with /3 — *■ 0, its varia- 
tion also tends to 0, so yF.p becomes an increasingly ac- 
ceptable semiclassical tunneling solution with an asymp- 
totic extremal action Se,f = f 71 " 2 - Returning to the 
original problem, I expect Se 

-1/2 



and y —> 

, K r'''?/F,^0i which is largely what I find. (For expo- 
nentially large |k| one expects R — > — In |«|; but the 
fitting function of Eq. (15) was chosen to obtain a simple 
and adequate fit, not to generate the correct asymptotic 
behavior.) 

The above results allow an easy calculation of the tun- 
neling action, and also of the bubble's field configuration 
for alH > and |af| > t by using the bounce solution 
with argument p = v 1 x 1 — t 2 . (Within the lightcone, the 
field typically oscillates about the true vacuum; for this 
post-tunneling evolution, the equations of motion must 
be solved anew with imaginary p.) 

One sample application of these results is |Q in the 
minimal model of supersymmetric gauge mediation, in 
which the large Higgs VEV hierarchy can generate a new 
global minimum in the potential of the scalar superpart- 
ners of the tau lepton. The potential in fact involves 
three fields, two of which are electrically charged and 
so should not acquire a VEV. Using the above methods, 
and restricting to a straight path connecting the false and 
true vacua, a significant range of that model's parame- 
ters is ruled out on the grounds that tunneling would be 
too fast. A small span of parameters remains for which 
a true multidimensional analysis would be needed to test 
whether the lifetime would be long enough. Note, how- 
ever, that no vacuum stability study can establish that a 
parameter value is definitely allowed: even if it produces 
a long-lived false vacuum, the evolution of the early uni- 
verse might put the initial state of <f> near the true vacuum 
and the model would be unacceptable. Vacuum stabil- 
ity can only rule out parameter values, and the present 
analysis does just that. 

Somewhat stronger bounds may result not only from 
tunneling via other paths through a multidimensional 
field space, but also from thermal fluctuations in the 
early, hot universe (for a related discussion, see Ref. jl| 
and references therein). However, such bounds are in 
general very model-dependent and, in the latter case, 
maybe subject to more uncertainties (see, e.g., Ref. [p(i[). 

The results of this work provide significant, robust and 
easy to determine bounds on a wide class of models in 
which the acceptable vacuum which we presumably in- 
habit can be destabilized by quantum tunneling. It is 
true that for such models the shape of the bounce, also 
determined in this work, is of only academic interest to 
us, since we will not survive the tunneling. The anal- 
ysis of this shape can be useful, however, in a different 
class of models: ones in which the true vacuum is the 
one we live in, while the phase transition occurred long 



ago. Then the profile of the bubble determines various 
properties of the early universe, such as how and where 
the latent heat was deposited, what was the spectrum of 
the particles produced, and what remnants were left be- 
hind. To properly study such a scenario may require the 
inclusion of gravity, the time-dependent evolution of the 
universe, and possibly nonzero temperature. The results 
of such studies will quantitatively be quite different from 
the present work, but hopefully much of the methodology 
introduced here will still be useful. 

Enlightening discussions with S. Coleman, I. Das- 
gupta, S. Kusenko, A. Linde and especially R. Rattazzi, 
who collaborated in an early stage of this work, are grate- 
fully acknowledged. 
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